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Abstract. 

Coherent coupling of two qubits mediated by a nonlinear resonator is studied. It 
is shown that the amount of entanglement accessible in the evolution depends both on 
the strength of nonlinearity in the Hamiltonian of the resonator and on the initial 
preparation of the system. The created entanglement survives in the presence of 
decoherence. 
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1. Introduction 

Developments in quantum information science rely critically on entanglement - quantum 
correlations between two or more physical systems e.g. between two qubits. There are 
various ways of creating entangled states of two qubits. The most natural is to place 
the qubits close to each other so that they can interact directly by local interactions i.e. 
a mutual inductance or capacitance. There are spectacular experiments performed for 
solid state qubits using this type of interaction [H [2[ [3l SI |5] . However the important 
stetp towards implementing quantum computation is a controllable coupling of qubits. 

The entanglement of distant qubits can be reached by swapping of entanglement 
[6l [TJ [8] or by some 'interaction-carrying' medium: a bus. In the following we 
discuss the latter case. One of natural candidates for such a bus are photons 
that are highly coherent. The problem of interaction of qubits with a quantized 
monochromatic electromagnetic field in a quantum cavity has been discussed in several 
papers P, [101 [HI [121 [E]- The field mediating the entanglement is usually described 
by a linear resonator [9l [TOl [Til [IH IIS]- However, it has been shown [11] that strong 
entanglement of qubits in linear coupling regime is possible only for specific initial states. 
In this paper we want to check wheather a nonlinear radiation improves this situation. 

Nonlinearity of the superconducting circuit resonator operating at microwave 
frequencies results in highly non-classical properties of the electromagnetic radiation 
generated by the resonator [16j. Such circuits coupled to the qubits has recently been 
intensively studied [17] since they do provide a natural measuring device [HI [19]. On 
the other hand, the problem of nonlinear oscillator mediating the interaction between 
qubits is still open for theoretical investigations. A pair of qubits coupled via a driven 
Duffing oscillator has been studied in [20] with an emphasis given on the (semi) classical 
limit of the oscillator. In some coupling schemes [211 [221 [23] an intermediate qubit with 
one or more Josephson junctions plays a role of a coupler. As the junction is strongly 
nonlinear, the coupler is treated as a nonlinear resonator. 

A nonlinear character of coupling has been also broadly exploited in quantum 
optics: there are many proposals where nonlinear coupling effects play an essential 
role in generation of entanglement of photons [2^ [251 [26l [27] and in construction of 
quantum gates [28j or teleportation protocols [29] . 

In this paper we assume two qubits strongly coupled to a resonator and discuss 
a coherent, non-local coupling between qubits via a linear and nonlinear resonator. 
The goal is to determine conditions which allow to strongly entangle the qubits with 
each other with the disentangled resonator. Thus the careful study of properties of the 
quantum bus is needed in order to fulfill this task. We investigate the infiuence of the 
strength of nonlinearity on coherent coupling of qubits and we show that it is highly 
nontrivial. On the basis of obtained results we propose a scheme of a tunable entangling 
gate. 

The results are applicable to entanglement through optical and microwave nonlinear 
radiation under the assumption that the qubit-field interaction is sufficiently weak to 
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be described in the rotating wave approximation [50] . 

In the first part of the paper we assume that all effects of dissipation and 
decoherence are negligible over studied time scales. We briefly discuss the effect of 
dissipation in the last section of the paper. To quantify the entanglement we calculate 
the negativity N [3l] which is a natural measure of entanglement. Our discussion 
is general enough to not depend on the specific architecture. However, a candidate for 
possible implementation is the fiux or charge qubit interacting with microwave radiation. 

In Section [2] we present the model of the system, in Sec. 13.11 and 13.21 we trace 
the evolution of the originally disentangled state with a single excitation. We consider 
two initial states which will be shown to result in qualitatively different behaviour of 
the qubit-qubit entanglement: nonlinearity-assisted entanglement gain or entanglement 
suppression. In Section H] we investigate the system with two and three excitations, in 
Sec. E] we verify the results for the negativity by presenting the qubit-qubit density 
matrix. Finally in Sec. [6] we shortly study the influence of decoherence on the qubits 
entanglement. Discussion and conclusions are given in Sec. [71 



2. The Model 



The system under considerations consists of two qubits Qi, i = 1,2 and a monochromatic 
resonator R. We assume that both qubits are coupled to the resonator simultaneously. 
To simplify the notation we apply units such that h = 1. In the absence of direct 
qubit-qubit interaction the Hamiltonian of such a system takes the form 

H = Hq^ + Hq^ + Hji + Hq^r + Hq^r, (1) 

the qubit Hamiltonian is 

Hq. = ycr., (2) 
The interaction terms Hq.r read 



where = + icTy, cr^ = {a^Y and cr^ are the Pauli matrices and 7 is the qubit-fleld 
coupling constant depending on the speciflc architecture. The resonator is described by 
a nonlinear bosonic oscillator of the Hamiltonian 

Hr = ujR (^a^a + ^) + Vr, (4) 

where the term Vr describes the nonlinearity. We shall discuss two different forms of 
the nonlinearity. The flrst is a polynomial one 

= a (a' + {a^y) (5) 

which is known to be related to the celebrated squeezed states [321 [33]. The second one, 
of central interest for our considerations, is the cosine-like nonlinearity 

= a cos(a + a^) (6) 
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The cosine term introduces nonlinearities of all orders and in that sense the cosine is 
'more nonlinear' than any polynomial. This type of nonlinearity is present e.g. in the 
Hamiltonian of the microwaves generated by a SQUID [H [231 El]- For such a circuit 
the parameter a can be changed by introducing into the resonator an adjustable weak 
link loop with two Josephson junctions. The Josephson tunneling energy controllable 
can by an external flux (p^, in this case a = 2Ej cos{7i(j)c/4>o) [Ml [35]. Varying 0c we 
can change a between 2Ej and 0. A tunable anharmonic LC circuit is also important 
as it introduces non-uniform level spacing reducing leakage to higher states [21] . In the 
limit of Vr = one arrives at the Hamiltonian of the Jaynes-Cummings type [36] for 
which exact solutions are known [TT] . 

The state vector of the system at t = is a tensor product of three constituent 
states 

\i^QiQ2Bit = 0)) = I^Q,) ® I^Q,) ® I^Pr) (7) 

This choice is certainly justified for weakly interacting systems. The unitary evolution 
of the initially factorisable state leads in general to the entangled tripartite state. By 
taking the trace over the R states one obtains the reduced density operator Pq^q^ which 
we then use to calculate the negativity [31] 

N{pQ,Q,)=max{0,-Y.Xi), (8) 

i 

where Aj are negative eigenvalues of the partially transposed [37j density matrix of the 
two qubits. is for separable states and reaches its maximal value N = 1/2 for 
maximally entangled states. 

We look for the conditions under which the qubits get very strongly entangled. 
3. Entanglement and nonlinearity 

In the following we study the influence of the nonlinear term in the bosonic Hamiltonian 
(jlj) on the quantitative and qualitative properties of entanglement of the qubit-qubit 
system. We consider below only the cosine nonlinearity as we have found that the 
influence of the nonlinearity ([5]) is weaker and not qualitatively different. We compare 
the results with those obtained for the linear single mode resonator, which despite 
its simplicity is a natural reference system. For numerical analysis we truncated the 
photonic space to M = 40 for which Tr{\ipQ^Q2R{t)){iljQ^Q2R{t)\ ~ 1). 

3.1. Entanglement gain 

The first of the considered initial states is 

l^QrQ.Bit = 0)) = \egO) (9) 

Here e and g stand for the excited and ground state of one qubit respectively. The state 
([9]) describes the first qubit in an excited state, the second one in the ground state and 
there are no photons inside the resonator. 
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Let us start with an analysis of the system in the absence of the nonhnear term i.e. 
Vr = 0. In this case one can perform exact analytical calculations of the wave function 
and trace the behavior of the system. The wave function of the system at the time t is 

m 

I^QiQ.i?W) = ^(l + cos(7t))|e^70)) (10) 

1 i ^ 

- - (1 - cos(7t)) l^eO) + -j=sm{^t)\ggl), 

where 7 = \/2'y. The time evolution goes in the subspace spanned by the states 
\egO),\geO),\ggl). For ^t = rni, n = 1,2,... for which the R subsystem could be 
decoupled from the qubits also the coefficient at If^eO) term becomes zero and the whole 
system remains disentangled. It is shown in Fig{T]as lines A {QQ — R entanglement) and 
B {QQ entanglement) at ut ~ 220. For other values of 7^ we are not able to obtain the 
entanglement of the qubits without the intrusion of the resonator states. This results 
in weak QQ entanglement with the limiting value Nqq^ ~ 0.1. The inclusion of the 




Figure 1. (color online) The QQ — R (line A) and QQ (line B) negativity as a function 
of dimensionless time Lut for the linear case and for the initial state |egO). Lines C and 
D depict the QQ-R and QQ negativity for the initial state \ggV)\ — lo,^ = O.Olw. 
Maximal qubit-qubit entanglement is possible when the resonator is not entangled with 
the qubits (lines C and D a.t uit — 110). 

nonlinear term {Vr 7^ 0) allows to go beyond this limit. It is shown in Figj2]that even 
for the relatively weak nonlinearity the negativity remains bounded by nothing but its 
natural limit i.e. N^^'' = 1/2. 

Increasing the strength of the nonlinearity (FigJSj) the qubits get permanently entangled 
for a relatively long time which increases with increasing a. For very strong nonlinearity 
{a = 2) the period is 25000ci;t and the maximal value of Nqq ?a 0.45 (not shown). Thus 
we see that for the non-symmetric initial state, when the qubits start from opposite 
states, the nonlinear resonator can lead, in contrast to the linear one, to the emergence 
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of strongly entangled QQ states remaining very weakly entangled with the resonator 
states. In other words we get the coherent quantum state transfer between the qubits 
through the quantum bus. 
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Figure 2. (color online) The unitary evolution of the qubit-qubit entanglement as a 
function of dimensionless time cot and the nonlinearity coefficient a. The initial state 
\egO). =uj, J = 0.01lj 
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Figure 3. (color online) The qubit-qubit negativity for large values of the nonlinearity 
coefficient a. The initial state \egO). fli = uj,'^ = O.OIlj. 



3.2. Entanglement suppression 

The situation looks different for the 'symmetric' initial state 



(11) 
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when the excitation is placed in the resonator and both qubits are in the ground state. 
The analytical calculations for the linear resonator (V/j = 0) result in the wave function 

m 

I^QiQ2flW) = [sin(7t)|ec/0) + sin(7t)|^eO)] 

+ cosm\ggl). (12) 

We see that at ^t = it/2 + mvr, where m is an integer, we get maximally entangled 
qubit-qubit (Bell) state 

1 



\Bi) 



V2 



{\e9) + \9e)). 



(13) 



It is important that the qubits can be in the maximally entangled state only, when their 
common state space is separable from the resonator space. It is shown in FigH] (lines C 
- the QQ-R negativity and D - the QQ negativity) at cjt = 110. As in the case of the 
initial state \egO), the qubit-qubit negativity Nqq reaches the value 0.1 (line D) when 
the QQ-R entanglement is maximal (FigH] cj/f: = 50 line C) . 

The inclusion of the nonlinear term spoils the symmetry (see Eq. (fT2!) ) and results in 
suppression of the qubit-qubit entanglement with the negativity limited by Nqq^ < 1/2 
depending on the amplitude a (see FigH]). 




1000 



Figure 4. (color online) The influence of nonlinearity strength a of the resonator on 
the entanglement of the qubits for the initial state (|lip . fli = uj,j = O.Olw. 



4. States with many excitations 

Finally, let us proceed to a consideration of the evolution of the originally disentangled 
states with many excitations. In Figs [5] and M we show the resulting qubit-qubit 
entanglement for the initial state with two and three excitations respectively 

li^Q^QMt = 0)) = \egl) (14) 
\^pQ,QMt = 0)) = \eg2). (15) 
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If the interaction between the qubits goes via a hnear resonator {a = 0) the negativity 
Nqq (and hence the QQ entanglement) is very small (solid line). In other words, 
one is not able to get rid of the resonator degrees of freedom without leakage of the 
information about the state of the system. The situation changes if the interaction goes 
via the nonlinear resonator (a 7^ lines in Figs [5] and [6]) and we obtain strong QQ 
entanglement. Thus in this case the influence of the nonlinearity of the "quantum bus" 
is favorable for the creation of the entanglement of the qubits. 
If we start from 

I^q^qMO)) = |eeO) (16) 

the interaction via the linear resonator does not entangle the qubits {Nqq = 0). It is in 
agreement with exact analytical calculations [H] . We have also checked it by calculating 
the concurrence (not shown) which is another entanglement measure. Surprisingly, 
switching "on" the nonlinearity does not improve the situation and the qubits remain 
disentangled. 

For the initial state 

1^(0)) = \gg2) (17) 

for the interaction via the linear resonator we get the QQ entanglement with Nqq^ ^ 
0.18. The presence of nonlinearity leads, in general, to a decrease of Nqq^ but for some 
specific a (a ~ 0.7) we have found a small increase of the negativity (not shown). 

Summarizing, the influence of nonlinearity is favorable for the initial states with 
one qubit in the excited state. Then the coupling via the quantum bus leads to strong 
entanglement and a coherent state transfer between the qubits. 
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Figure 5. (color online) The influence of the nonlinearity strength a on the qubit- 
qubit entanglement for the initial state \egl), Qt ~ uj,^ — O.Olw. 
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Figure 6. (color online) The influence of the nonlinearity strength a on the qubit- 
qubit entanglement for the initial state \eg2), D,i = lo,^ = O.Olw. 
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Figure 7. (color online) The real part of the qubit-qubit density matrix for the case 
of the linear resonator and the initial state ([9]). The matrix is taken at tot = 111 for 
which Nqq = 0.104 (solid line in Figl2]). The imaginary elements of the matrix are of 
the order of lO^^**. 



5. Qubit-qubit density matrix 

One of the best demonstrations of quantum entanglement is the quantum state 
tomography which yields a density matrix of coherently coupled qubits [38]. Such 
experiments for qubits interacting via a resonant cavity are in progress [13]. As an 
example we present results of our calculations of the qubit-qubit density matrix for the 
case discussed in Sec J3.1l 

In FigJT] we show the real part of the density matrix for the qubits interacting via 
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Figure 8. (color online) The real part of the qubit-qubit density matrix for the case 
of the nonlinear resonator (a = 0.0035) and the initial state The values are taken 
at ujt = 435, Nqq = 0.492 (dashed line in FigEJ . 





Figure 9. (color online) The imaginary part of the qubit-qubit density matrix for the 
case of the nonlinear resonator {a ~ 0.0035) and the initial state ([Sj). The values are 
taken at ujt = 435, Nqq = 0.492 (see dashed line in FigH]). 



the linear resonator. The matrix elements are taken at out = 111 for which Nqq = 0.104 
in Fig 121 Looking at the state vector ( |T0|) . this situation corresponds to equal probability 
of finding the excitation in the resonator or in one of the qubits. In FigjH] and Figj9] we 
present the real and imaginary part of the QQ density matrix for a = 0.0035. We see 
that in this case the non-diagonal matrix elements being the hallmark of entanglement 
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have much larger values as when the coupling goes via the linear resonator. The values 
correspond to dashed line in Fig|2]at out = 435, Nqq = 0.492. 

6. The effect of decoherence 

Design and construction of quantum devices is always affected by the influence of 
environment. In the following we apply the commonly used Markovian approximation 
[30] and model the reduced dynamics of the QQR system in terms of master equation 
generating complete positive dynamics [39]. We take into account the decoherence of 
both qubits and the resonator. Following [9l [8] we assume that the effect of environment 
can be included in terms of two independent Lindblad terms: 

p(t) = [Lj, - hn]pit) (18) 



where the 'conservative part' is given by 

LHi-) = -t[H, ■] 
whereas the 'Lindblad dissipator' 



i=l 



Lr{ 



A^A{-) + {■)A^A-2A{-)A^ 
S,tS,(-) + (-)SjS,-2S,(-)Sl 



(19) 



(20) 

(21) 
(22) 



is expressed in terms of creation and annihilation operators 'weighted' by suitable 



lifetimes A = a/ \/Tr and S j = cr / ^ 



where T^^ = 5 ■ 10 and Tr 



10" 



are the resonator and the qubits decoherence times respectively. It is shown in FigJTOl 
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Figure 10. (color online) The influence of decoherence on the entanglement of the 
qubits for the initial state ©. f]^ = cj,7 = 0.01w,a = 0.0035, Tr, 5 • 10"^ Tq^ = 
10-S. 



that the effect of such a local damping is purely quahtative and the constructive role 
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of of the nonlinear term in the oscillator is not obscured provided that the relaxation 
times are sufficiently large. 

7. Discussion and conclusions 

In this paper we have considered entanglement generation in the system consisting of 
two qubits coupled by a linear or nonlinear resonator. The goal was to test the role of 
nonlinearity and to find conditions, under which the qubits became strongly entangled. 
We have considered cosine-like nonlinearity where the discussed effects become most 
visible. We have found that the infiuence of the nonlinearity of the resonator, having 
its own dynamics, depends significantly on the initial state of the investigated system. 
We have shown two qualitatively different behaviours. The first occurs in the system 
with non-symmetric initial states. The presence of nonlinearity results in strong 
enhancement of the entanglement almost up to its maximal value. For sufficiently large 
amplitude of the nonlinear term, the entanglement becomes (quasi) permanent and does 
not vanish for a long time. This constructive role of nonlinearity is certainly desired for 
the applications. For the symmetric initial states the role of the nonlinearity is no more 
constructive but rather results in suppression of entanglement. 

On the basis of these findings, we propose the following tunable coupling scheme: 
for the non-symmetric initial state the nonlinear term should be switched on, for the 
symmetric one this term should be Oj^^and the interaction goes via the linear resonator. 
Then in both cases we obtain the desired strong entanglement of qubits and the system 
acts as a quantum entangling gate. To our knowledge, such a scheme has not been 
proposed yet. 

The main advantage of this gate is that the strong entanglement of qubits can 
be reached for many initial states, even for the states with many excitations (with 
some exceptions). We have also shown that the coherent coupling of qubits survives 
in the presence of dissipation assuming decoherence times in agreement with recent 
experiments [T2l [T3] . 
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